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The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect 
the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface 
detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies 
and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly 
calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per 
area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional 
function that takes into account signal asymmetries due to interference between the geomagnetic and charge- 
excess emission components. The spatial integral over the signal distribution gives a direct measurement of the 
energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 
15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This 
radiation energy - corrected for geometrical effects - is used as a cosmic-ray energy estimator. Performing an 
absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator 
scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution 
of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with 
at least five radio stations with signal. 
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I. INTRODUCTION 

Cosmic rays in the ultra-high energy regime are detected 
through giant particle showers developing in Earth’s atmo¬ 
sphere. Various detection systems are used to measure the 
calorimetric energy of the shower. Well established are tech¬ 
niques using telescopes that observe directly the shower de¬ 
velopment through fluorescence light emitted by molecules 
excited by the shower particles, and/or detectors positioned 
on the surface of Earth that measure the particles at one stage 
of the air-shower development (e.g., [1, 2]). 

The observation of radio signals emitted by the shower par¬ 
ticles using broadband megahertz (MHz) antenna stations has 
also been explored as a detection method to obtain comple¬ 
mentary information on the air-shower development, and has 
become an active field of research in recent years [3-5], The 
properties of the primary cosmic rays have been studied in 
this way including their arrival direction, energy, and com¬ 
position. Directional information can be obtained from the 
arrival times in several radio stations [3, 4, 6, 7]. To obtain 
information about the energy, calibrated detectors for cosmic- 
ray showers co-located with the radio stations are used [7-9]. 
Composition information has also been derived by relying on 
simulations of radio emission [9, 10]. 

One of the interesting characteristics of the radio-emission 
signal is the strong polarization of the electric field arriving at 
the antennas. Two components have been identified originat¬ 
ing from different emission processes. The dominant one is 
perpendicular to Earth’s magnetic field and is denoted as ge¬ 
omagnetic emission [3, 6, 11], The second component is po¬ 
larized radially with respect to the axis of the air shower and 
results from the negative charge excess in the shower front 
[12-14], Its relative strength with respect to the geomagnetic 
emission is on average 14% at the Auger site for an air shower 
arriving perpendicularly to the geomagnetic field [15]. 

As a consequence of the superposition of the two emission 
mechanisms, the lateral distribution function (LDF) of the 
electric-field strength has been found to have a radial asym¬ 
metry [16-19]. The two-dimensional shape of the LDF is best 
understood in a coordinate system with one axis perpendicular 
to the shower direction v and Earth’s magnetic field B (along 
the Lorentz force ~ v x B), and another along the perpen¬ 
dicular axis v x (v x B). In this coordinate system the LDF 
exhibits a peanut-like shape. 

So far, all radio experiments have used experiment-specific 
quantities to reconstruct the cosmic-ray energy, such as the 
radio signal strength at a characteristic lateral distance from 
the shower axis. While this method has long been known 
to provide a good precision [20], it has the disadvantage that 
the corresponding energy estimators cannot be directly com¬ 
pared across different experiments. The main reason for this is 
that the shape of the lateral signal distribution changes signif¬ 
icantly with observation altitude. The optimal characteristic 
distance varies with observation height and even at the same 
characteristic distance the radio signal strengths are signifi¬ 
cantly different [21]. Hence, a comparison between different 
experiments cannot be performed directly. 


In this contribution we introduce a general approach with 
a direct physical interpretation. At each observer position we 
calculate the energy deposit per area of the cosmic-ray radio 
pulse and by integrating the two-dimensional lateral distribu¬ 
tion function over the area we obtain the total amount of en¬ 
ergy that is transferred from the primary comic ray into radio 
emission during the air-shower development. This approach 
is independent of the shape of the signal distribution because 
energy, i.e., the integral over the signal distribution, is con¬ 
served. 

In this analysis we present the relation between the cosmic- 
ray energy and the total energy emitted by the air shower as a 
radio pulse, for primaries of energy in the EeV (= 10 18 eV) 
range. To obtain this relation we use radio stations of the 
Auger Engineering Radio Array located within the Pierre 
Auger Observatory in Argentina. The antennas are of the log¬ 
arithmic periodic dipole antenna (LPDA) type and have been 
thoroughly studied and calibrated [22]. We take advantage of 
the possibility to cross-calibrate these measurements with the 
well-understood data of the Observatory and with recent de¬ 
velopments in understanding the radio-emission mechanisms, 
with their corresponding polarization patterns of the electric 
field and the particular lateral distribution of the total field 
strength [19, 23]. 

This paper is structured as follows. We begin with the ex¬ 
perimental setup of the antenna array and the surface detec¬ 
tor, then proceed to the data selection and event reconstruction 
procedure. After that we describe the calibration that uses a 
likelihood procedure, and we discuss experimental uncertain¬ 
ties. Finally, we present the energy measurement of the AERA 
radio detector, its resolution, the correlation of the radiation 
energy with the shower energy and we address the systematic 
uncertainty of the radio energy as an energy estimator. 

A summary of the main results presented here and its im¬ 
plications on the energy measurement of cosmic rays can be 
found in an accompanying publication [21]. 


H. DETECTION SYSTEMS 

The Pierre Auger Observatory is a hybrid detector for 
cosmic rays, based on two complementary detection sys¬ 
tems. The surface detector (SD) array consists of 1660 water- 
Cherenkov detectors distributed over an area of 3000 km 2 . Its 
stations have a spacing of 1.5 km, optimized to reach full effi¬ 
ciency for cosmic-ray energies above 3 EeV [1], The fluores¬ 
cence detector (FD) consists of 27 telescopes grouped at four 
locations around the area covered by the SD stations. With the 
FD, UV light is observed originating from the fluorescence 
emission of molecules excited by the cosmic-ray-induced air 
shower. The hybrid design of the Pierre Auger Observatory 
allows for an accurate energy calibration of the SD using the 
direct energy measurement of the FD. The amount of fluores¬ 
cence light is proportional to the deposited energy and thus 
yields an accurate measurement of the energy of the primary 
particle. 

The radio detector (RD) stations of AERA are located in an 
area of denser detector spacing of the SD array. This region. 
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with SD station spacing of 0.75 km, allows the detection of 
cosmic-ray energies down to about 0.1 EeV. 

The first deployment stage of AERA consists of 24 antenna 
stations with a spacing of 144 m. Every station is equipped 
with two logarithmic-periodic dipole antennas [ 22 ] integrated 
in one mechanical structure. The two antennas are oriented 
into the east-west and north-south directions relative to mag¬ 
netic north. The corresponding analog and digital electronics 
are tuned to the frequency range of 30 to 80 MHz [24]. After 
filtering and amplification, the signal is digitized at 180 MSa/s 
or 200MSa/s depending on the hardware type [25]. The sta¬ 
tions are equipped with solar cells and a battery to ensure an 
autonomous power supply. Furthermore, all stations are con¬ 
nected via an optical fiber-network to the data-acquisition sys¬ 
tem (DAQ). 

The system runs in two different modes, depending on the 
type of digitizing hardware. A self-trigger algorithm runs on 
the voltage trace itself, which identifies pulses based on char¬ 
acteristics described in [26], and consequently creates a trig¬ 
ger. The triggers of multiple RD stations are checked for co¬ 
incidences in a short time window of 1 qs - compatible with 
the passage of an air shower - at the DAQ level. A readout is 
requested once coincidences between at least three radio sta¬ 
tions are found. Alternatively, stations are triggered using an 
external trigger. Here, the DAQ receives a trigger from the 
Observatory’s central data-acquisition system (CDAS) once 
an air-shower candidate has been registered with the SD or 
FD. This trigger initiates the readout of all the stations, which 
are equipped with a ring buffer. The buffer has a size of 4 
GBytes and can store the traces of the two channels for about 
7.4 s which is sufficient to hold the data for the time needed to 
receive the trigger by the CDAS. 

III. DATA SELECTION AND EVENT RECONSTRUCTION 

In this work we are using RD and SD data recorded between 
April 2011 and March 2013 when AERA was operating in 
its first commissioning phase. The data are stored as events, 
which refer to all relevant information that has been read out 
following a trigger. For this analysis, both self-triggered and 
externally triggered events are used. 

A. Preselection of cosmic-ray candidates 

In the case of the self-triggered events, a preselection is per¬ 
formed offline by searching for coincidences with the surface 
detector events. A radio event has to agree in time and location 
with an SD event to be considered as cosmic-ray candidate. 
The radio-trigger time and the time when the air shower core 
hits the ground have to agree within ±20 qs. Such a conserva¬ 
tive coincidence window also accounts for horizontal events, 
for which the time difference is expected to be larger. 

For both trigger types, only events with a clear radio pulse 
in at least three stations are considered, to allow for a recon¬ 
struction of the incoming direction of the signal. For exter¬ 
nally triggered events the requirement is a signal-to-noise ra¬ 


tio (SNR) greater than ten. Here the SNR is defined as the 
maximum of the Hilbert envelope-squared [27] divided by the 
noise variance. For self-triggered events the signal threshold 
is dynamically adjusted to the noise level to keep the trigger 
rate at a constant level of 100 Hz. We require that the recon- 
structed incoming directions from the radio and the surface 
detectors agree within 20 ° to be accepted as a cosmic-ray can¬ 
didate. The 20° cut does not reflect the angular resolution of 
the SD nor that of the radio detector. This preselection cut re¬ 
tains the maximum number of cosmic-ray signals and signif¬ 
icantly reduces the number of random (anthropogenic) noise 
pulses, which originate mainly from the horizon. 

In addition, we apply quality cuts on the data of the surface 
detector [28]. The most important cuts are that the core posi¬ 
tion is closest to an active station and surrounded by a hexagon 
of active stations and that the zenith angle of the incoming di¬ 
rection be less than 55°. A total of 181 cosmic-ray candidates 
with energies above 10 17 eV remain. 

As an engineering array, AERA was subject to several 
changes in software and hardware which significantly limited 
the uptime. In future, we expect a larger rate of cosmic rays 
due to the stabilized operation of the detector. 

B. Reconstruction of radio data 

We use the software framework Off line [29] of the Pierre 
Auger Collaboration to process the measured raw data. First, 
the air shower is reconstructed using the surface detector in¬ 
formation [30]. Second, the reconstruction using the radio de¬ 
tector data is performed [31]. Narrowband noise sources are 
filtered out using a radio-frequency interference suppression 
in the time domain. Sine waves with the frequency of noise 
sources are fitted to the measured voltage trace and subtracted. 

We correct for the influence of the analog signal chain using 
the absolute calibration of the AERA station and reconstruct 
a three-dimensional electric field by using the direction of the 
shower and applying the simulated antenna response [ 22 ]. 

An example of a reconstructed electric-field trace E(t) is 
shown in Fig. 1. The energy fluence /, i.e., the energy deposit 
per unit area, of the incoming electromagnetic radio pulse at 
each radio station is determined by calculating the time in¬ 
tegral over the absolute value of the Poynting vector. This 
is achieved by squaring the magnitude of the electric-field 
trace and summing over a time window of 200 ns ([( 1 ,( 2 ]) 
around the pulse maximum which has been determined from 
the Hilbert envelope of the trace (cf. Fig. 1). The contribution 
of background noise (determined in the noise window [( 3 , £ 4 ]) 
is subtracted under the assumption that the main contribution 
is white noise. The energy fluence / is given by 

/ = e 0 c (AtJ2\E{ti)\ 2 - _ j- \Eiti)\ 2 j , (1) 

\ tl £3 / 

where Sy is the vacuum permittivity, c is the speed of light in 
vacuum and At is the size of one time bin. This quantity is 
used throughout the whole analysis and will be given in units 
of eV/m 2 . To approximate the uncertainty the noise level as 
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FIG. 1. Reconstructed electric-field trace of one of the measured 
cosmic-ray radio events. An upsampling by a factor of five was ap¬ 
plied. The shown Hilbert envelope (dashed line) is the square root of 
the quadratic sum of the Hilbert envelopes of the three polarization 
components. 

described above is used. As the radio detector effects have 
been corrected for, the reconstructed energy fluence can be 
directly compared to air-shower simulations. 

We also calculate the direction of the electric-field vec¬ 
tor, i.e., the polarization direction of the signal. In the full 
width half maximum (FWHM) interval around the pulse max¬ 
imum of the Hilbert envelope we observe that the recon¬ 
structed electric-field vectors are aligned approximately along 
the same direction for every time bin. To accurately determine 
the mean direction of the electric-field vector, we average over 
all vectors in the FWHM interval of the Hilbert envelope (cf. 
Fig. 1). 


C. Selection of radio signals induced by cosmic rays 

Given the amount of pulsed background noise at the AERA 
site, the preselected events are likely to contain non cosmic- 
ray signals that mimic cosmic-ray pulses. There are two sce¬ 
narios possible: Signals in one or more stations are not caused 
by the air shower or an event contains only noise pulses that 
by chance led to a reconstructed incoming direction similar to 
that of the SD. 

In order to reject background signals, we take advantage of 
the expected polarization of the radio signal. The polariza¬ 
tion of the radio pulse is only used for this purpose and not 
considered for the energy estimation. In the frequency range 
of AERA (30 to 80 MHz) the dominant emission process is 
the geomagnetic emission [11, 15]. Here, a linear polarization 
of the electric field is expected to be in the direction of the 
Lorentz force (given by e geo ) that acts on the charged particles 
while they traverse the magnetic field of Earth. The polariza¬ 
tion is altered by an additional emission which is linearly po¬ 
larized radially towards the shower axis (given by ecE), and is 
referred to as the charge-excess emission process [15, 32-34], 

The expected direction of the electric-field vector is there¬ 
fore calculated from the geomagnetic and the charge-excess 
contributions 

E ex p oc sin a e geo + a e c e , (2) 



FIG. 2. Polarization map of a single event. The axis coordinates 
are in the shower plane where the x-axis corresponds to the direction 
of the Lorentz force (v x B) and the y-axis perpendicular to that 
and to the shower axis (v x (v x B)). The SD shower core is at the 
coordinate origin. The measured polarizations are shown as the black 
arrows. The gray arrows are the model expectations, and the red and 
blue arrows are the geomagnetic and the charge-excess components, 
respectively. The definition of /? is described in the text. The air- 
shower properties of this event are: Energy of 0.9 EeV and arriving 
from a zenith angle of 36° and from 27° south of west. For the 
emission model of Eq. (2), the optimal value of the relative charge- 
excess strength is a = 0.18. 

where a is the angle between the shower axis and magnetic 
field of Earth, and a is the average relative charge-excess 
strength that has been measured to be 0.14 ± 0.02 at AERA 
[15], In this approach, the direction of the geomagnetic con¬ 
tribution depends only on the incoming direction of the air 
shower whereas the charge-excess contribution depends in ad¬ 
dition on the position of the radio station relative to the shower 
axis. 

In Fig. 2, all stations with signal of a cosmic-ray candidate 
are shown, and the measured polarization is compared with 
the expectations of the two radio-emission mechanisms. The 
overall agreement between measured and expected field po¬ 
larizations is quantified using the angular difference 

Pi 4^(E/ meas? 2, -E'expji) (3) 

at each station i. For each event, the average deviation p of 
the individual deviations pi of the stations with signal is cal¬ 
culated and will be used as criterion for a quality cut. Relevant 
uncertainties are taken into account as follows: 

• The relative strength a of the charge excess can vary 
from event to event due to shower-to-shower fluctua¬ 
tions, and additional dependencies on the geometry of 
the air shower [35]. Therefore, for each possible values 
of a between 0 and 0.5 the average deviation P is cal¬ 
culated and only the smallest value of P is considered. 
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• The uncertainty of the SD shower core position is taken 
into account by variation of the core within its estimated 
uncertainties. In our data set the uncertainty varies 
between 10 m and 80 m depending on the energy and 
zenith angle. For each trial of the core position f3 is 
calculated. Again, only the smallest value of B is con¬ 
sidered. 

• Interference of the cosmic-ray radio signal with noise 
pulses can alter the polarization. Simulation studies 
showed that for a single radio station the uncertainty 
in B d ue to noise is below 8° at detection threshold, and 
decreases to 1° at high signal-to-noise ratios. To obtain 
the average value of /3 for all radio stations in the event 
we compute a weighted mean with weights Wi = 1 /cr 2 . 
with o'/ 3 i being the expected uncertainty from the simu¬ 
lation. 

We impose a limit on the average deviation /3 of the polar¬ 
ization direction. This maximum deviation is fixed at a value 
of 3°. This value is slightly above the combination of the fol¬ 
lowing effects. 

The incoming direction of an air shower reconstructed with 
the surface detector has an uncertainty between 1.3° and 0.7° 
depending on the cosmic-ray energy and the zenith angle [30]. 
Hence, the expected direction of the electric-field vector will 
have the same uncertainty. All antennas are aligned to the 
magnetic north (or perpendicularly to the magnetic north in 
case of the other polarization direction) with a precision of 
better than 1° [36]. All antennas are uniformly constructed 
and the two antennas of a radio station are identical. Asymme¬ 
tries in the ground conditions have only negligible influence 
as the LPDA antenna is mostly insensitive towards the ground. 
A measurement at AERA has shown that the responses of all 
antennas differ by less than 0.3% [37], 

A difference in the amplification of the signal chain of the 
north-south and east-west polarized antenna will influence the 
polarization measurement. From an individual measurement 
of the signal chain of all antennas the uncertainty is estimated 
to be 2.5% which results in a polarization uncertainty below 
0.7°. 

In addition, we neglect the dependence of the relative 
strength a of the charge excess on the distance between ob¬ 
server position and shower axis [35, 38], For a single station 
this effect is relevant. However, in our approach we only use 
the average deviation of all stations with signal also taking 
into account the uncertainty in the core position. Therefore 
the distance dependence will mostly average out. We estimate 
that the remaining additional scatter is 1.5°. 

We account for individual radio stations being contami¬ 
nated with substantial noise signals by iterating through all 
configurations with only one and then more stations removed, 
down to the minimum of three stations. An event where the 
weighted average deviation j3 is greater than 3° for all station 
combinations is rejected. If fi is less than 3° for any station 
combination and the fraction of selected stations is larger than 
50% of the total number of stations with signal, the event can¬ 
didate is considered a cosmic-ray event and only the stations 
from this particular combination are used. After this cut 136 


events remain. The number of excluded single stations and 
complete events is compatible with the measured rate of noise 
pulses. 

Most of the events recorded during thunderstorm condi¬ 
tions appear to be rejected by this selection procedure as the 
strong atmospheric electric fields of a thunderstorm influence 
the radio emission and alter the polarization of the radio sig¬ 
nals [39, 40]. For two thirds of the events, a measurement of 
the atmospheric electric field is available. These events are 
checked for thunderstorm conditions using an algorithm de¬ 
scribed in [41]. Based on this check, two additional events 
were rejected. All cuts are summarized in Table I. 

TABLE I. Overview of selection cuts and the number of events sur¬ 
viving these cuts. Preselection means: Ecr >0.1 EeV, standard SD 
quality cuts, > 3 radio stations with signal, SD and RD reconstructed 
incoming directions agree within 20°. See text for details. 


cut 

number of events af¬ 
ter cut 

preselection (Sec. Ill A) 

181 

polarization cut (/3 <3°, Sec. Ill C) 

136 

no thunderstorm conditions (Sec. Ill C) 

134 

LDF fit converged (cr < 300 m, Sec. IV) 

126 

> 5 stations with signal 


(only high-quality data set. Sec. V) 

47 


D. Uncertainties on the energy fluence in a single radio station 

In addition to the uncertainties on the amplification of the 
signal chain of 2.5% discussed above, no further uncertainties 
are expected that would result in a different response of sta¬ 
tions within one event. To first order, the frequency content 
and the incoming direction of the radio pulse are similar at all 
observer positions. Therefore, an uncertainty of the antenna- 
response pattern has a negligible influence as it is evaluated 
for the same direction at all stations. Possible different ground 
conditions at different station positions that result in a differ¬ 
ent reflectivity of the soil are negligible due to the insensitivity 
of the antenna towards the ground. The 2.5% amplification 
uncertainty results in 5% uncertainty on the energy fluence 
/, as / scales quadratically with the electric-field amplitude. 
This uncertainty is added in quadrature to the signal uncer¬ 
tainty resulting from noise. 


IV. ENERGY ESTIMATOR 

To obtain an absolute energy estimator from the signals 
at the different distances to the shower axis (energy fluence 
/ in units of eV/m 2 ) a LDF is used which takes into ac¬ 
count the signal asymmetries due to constructive and destruc¬ 
tive interference between the geomagnetic and charge-excess 
components, as well as Cherenkov time-compression effects 
[19]. This LDF describes the main features seen in simulated 
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FIG. 3. Lateral signal distribution of a single cosmic-ray event. The air-shower properties of this event are: Energy of 0.75 EeV and arriving 
at a zenith angle of 37° and from 44° west of south. Left: The energy fluence in the shower plane. The measurements are indicated as 
circles where the color shows the energy fluence. Gray squares are stations with signal below threshold and the red cross marks a station 
that is rejected due to a mismatch in the signal polarization. The background map shows the LDF parametrization. The coordinate origin 
is the reconstructed core position of the radio LDF fit. Note the lack of color contrast between the infill color of the data points and the 
background. This is indicative of the agreement between the data and the model. Right: Representation of the same data as a function of 
distance from the shower axis. The colored and black squares are the measured energy fluences and gray squares are the stations with signal 
below threshold. For the three data points with the highest energy fluence, the one-dimensional projection of the two-dimensional LDF onto 
lines connecting the radio-core position with the corresponding radio detector positions is illustrated with colored lines. This demonstrates the 
azimuthal asymmetry and complexity of the two-dimensional lateral distribution function. The inset figure shows the azimuthal direction of 
the three LDF projections. The distribution of the residuals (data versus fit) is shown as well. 


and measured cosmic-ray radio events. The LDF function is 
parametrized as 


f(r)=A 


exp 


( r + e vxB r core) 


— Cq exp 


( r T @2 e fixB r core) 

{C^e Ci a ) 2 


(4) 


All coordinates are in the shower plane, r denotes the sta¬ 
tion position. The four fit parameters are the amplitude A, 
the slope parameter a and the particle core position f core . In 
case of low station multiplicity, the particle core position is 
taken from the SD reconstruction, which enables us to also 
use events with only three or four stations with signal. Co - C4 
are constants that are estimated from CoREAS Monte Carlo 
simulations [23] and can be found in Appendix A. Co — C2 
are zenith-angle dependent. The LDF is fitted to the data us¬ 
ing a chi-square minimization. An example of one air shower 
within our data set is shown in Fig. 3. 

Some events do not contain sufficient information to fit the 
LDF, such as when only three stations with signal are present 
that have roughly the same signal strength. This results in 
an unphysically broad LDF. To reject these events we impose 
the quality cut a < 300 m (Table I). An analysis of air-shower 
simulations for the AERA geometry showed that the a param¬ 
eter of the LDF is never larger than 300 m. 

In the following, only the 126 events that pass the quality 
cuts are considered and will be referred to as the full data set. 
To derive the accuracy of the energy estimation method, the 


data set will be further divided in a high-quality data set con¬ 
taining only events with at least five stations with signal, i.e., 
events where the core position can be reconstructed in the ra¬ 
dio LDF tit. 


A. Definition of the energy estimator 

The spatial integral of the lateral distribution function gives 
the amount of energy that is transferred from the primary cos¬ 
mic ray into radio emission in the AERA frequency band dur¬ 
ing the air-shower development, and will be given in units 
of eV. We define the energy estimator ,S' ra( | 1() as this radiation 
energy divided by sin 2 a to account for different emission 
strengths at different angles between shower axis and mag¬ 
netic field, see Eq. (2), 

Sradio = . \ [ f{r)d 2 r 

sin a J 

K 2 (5) 

— An / 2 n ri'2 2C 4 <x\ 

— .9 \ a ^0^3 e ) 1 

sin a 

where Si 2 denotes the shower plane. The positive a 1 term 
dominates by far over the negative second term resulting in a 
positive value of Sj-adio- The sin 2 a correction only holds if 
the geomagnetic emission is the dominant contribution which 
is the case for <x > 10° at AERA. Due to the reduced emission 
strength the number of detections for arrival directions within 
10° of the geomagnetic field axis is suppressed. The angular 
distribution of the events is shown in Fig. 4. 





























9 


90 ° 



FIG. 4. Skymap of the 126 selected events. Green filled circles de¬ 
note air showers with at least five stations with signal and open cir¬ 
cles denote air showers with less than five stations with signal. The 
red star denotes the direction of the magnetic-field axis at AERA. All 
measured events are at least 20° away from the magnetic-field axis. 
Therefore, the geomagnetic emission gives the dominant contribu¬ 
tion to the radiation energy for all events. 

B. Event-by-event uncertainties of the energy estimator 

The following uncertainties are relevant for the energy es¬ 
timator due to event-by-event fluctuations and summarized in 
Table II: 

• The gains of the low-noise amplifiers and filter ampli¬ 
fiers exhibit a temperature dependence. The effect has 
been measured and amounts to — 42mdB/K. Each air 
shower is measured under specific environmental con¬ 
ditions. In particular this implies that we have a ran¬ 
dom distribution of ambient temperatures which exhibit 
a Gaussian distribution with a standard deviation of 8.3° 
C. This corresponds to a fluctuation of the gain of 4%. 

• An uncertainty of the simulated antenna response that 
depends on the incoming direction of the radio signal 
will lead to an event-by-event uncertainty as each event 
has a different incoming direction. The effect is deter¬ 
mined to be 5% by comparison of the simulated antenna 
response with a measurement at AERA [22], 

• The reconstmcted direction of the air shower obtained 
with the SD has an uncertainty of less than 1.3°. This 
has negligible influence on the antenna response pat¬ 
tern, since it can be considered uniform over such a 
small change of angle. 

As the different uncertainties are independent, the total uncer¬ 
tainty of the electric-field amplitude is a/4% 2 + 5% 2 ss 6.4% 


and therefore 12.8% on ,S' ra( i m . The uncertainty of a can be 
neglected. The fit uncertainties of A and a including their cor¬ 
relation are propagated into S' ra dio using Gaussian error prop¬ 
agation. In the case of events with less than five stations with 
signal, the core position of the surface detector reconstruction 
is used and its uncertainty is propagated into the fit uncertainty 
of Sradio- This fit uncertainty is added in quadrature to the 
statistical uncertainty of 12.8% of the energy estimator. The 
average fit uncertainty of 5 ra dio is 46%. For events with at 
least five stations with signal the average uncertainty reduces 
to 24%. 

C. Absolute scale uncertainties of the energy estimator 

The dominant systematic uncertainties of the reconstructed 
electric-field amplitudes are the calibration of the analog sig¬ 
nal chain and the antenna response pattern. The analog signal 
chain consists of the low-noise amplifier, the filter amplifier 
and all cables between the antenna and the analog-to-digital 
converter. The analog signal chain has been measured for each 
channel of each radio station separately in the field and dif¬ 
ferences are corrected for. The systematic uncertainty of the 
analog chain amounts to 6 %. 

The simulated antenna response pattern has been confirmed 
by measurements at an overall level of 4%. The systematic un¬ 
certainty of the measurement is 12.5% in the vector effective 
length [42]. Conservatively, the systematic uncertainty of the 
antenna-response pattern is therefore estimated as 12.5%. 

Systematic uncertainties introduced by the usage of the 
two-dimensional signal distribution function of Eq. (4) are 
negligible. Detailed comparisons of the shape of the radio sig¬ 
nal distribution measured with LOFAR with the predictions 
from CoREAS show no indication of any systematic discrep¬ 
ancy [43]. We determined the influence of the 2D-LDF model 
on the radiation energy in a representative CoREAS Monte 
Carlo data set for the AERA detector and found a systematic 
effect of less than 5%. 

Combining all uncertainties in quadrature, the systematic 
uncertainty of the electric-field amplitude is 14%. The radio¬ 
energy fluence and the energy estimator scale with the am¬ 
plitude squared. Therefore, the systematic uncertainty of the 
absolute scale of the radiation energy is 28%. We note that, as 
the cosmic-ray energy is proportional to the square root of the 
radiation energy (see next section), the systematic uncertainty 
of a radio cosmic-ray energy scale would remain at 14%. 

V. ENERGY CALIBRATION 

The radio-energy estimator S^adio is shown as a function of 
the cosmic-ray energy E’ch measured with the surface detec¬ 
tor in Fig. 5 top. A clear correlation is observed. For the cali¬ 
bration function we follow the same method as used for the 
calibration of surface detector events with fluorescence de¬ 
tector events of the Pierre Auger Observatory [44—46]. The 
calibration function 

S radlo = A x 10 7 eV (Ebr/lO 18 eV) B 


( 6 ) 
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TABLE II. Overview of uncertainties of the electric-field amplitude and the energy estimator S ra dio- “©” denotes a quadratic sum. The 
average fit uncertainty of S ra dio is 46%, and 24% for the high-quality subset of events with at least five stations with signal. 


source of uncertainty 

a \E\ 

^"■S’radio 

event-by-event 

temperature dependence 

4% 

8% 

angular dependence of antenna response pattern 

5% 

10% 

reconstructed direction 

negligible 

negligible 

LDF fit uncertainty 

- 

error propagation of fit parameters 

total event-by-event uncertainty 

6.4% 

12.8% © fit uncertainty 

absolute scale 

absolute scale of antenna response pattern 

12.5% 

25% 

analog signal chain 

6% 

12% 

LDF model 

<2.5% 

<5% 

total absolute scale uncertainty 

14% 

28% 


is obtained by maximizing a likelihood function that takes 
into account all measurement uncertainties, detector efficien¬ 
cies and the steeply falling energy spectrum (the functional 
form of the likelihood function can be found in appendix 
B). The result of the calibration fit is A = 1.58 ± 0.07 and 
B = 1.98 ± 0.04. The correlation between A and B is 35%. 
The resulting slope is quite compatible with an exponent of 
B = 2 implying that the energy deposited in radio emission 
increases quadratically with the cosmic-ray energy. If B is 
fixed to 2 the fit result is A = 1.59 ± 0.06. We can infer 
from Eq. (6) that, for a 1 EeV air shower perpendicular to the 
magnetic field axis, 15.8 MeV is deposited on average in radio 
emission in the frequency range of 30 to 80 MHz. 

The lower left panel of Fig. 5 shows the scatter around 
the calibration curve for all air showers in our data set. This 
amounts to 29%. We also tested a high-quality data set con¬ 
taining only air showers with at least five stations with sig¬ 
nal, where a determination of the core position in the radio 
LDF fit is possible. These air showers are marked by green 
filled circles in Fig. 5. The fit of the calibration curve gives a 
compatible result (A = 1.60 ± 0.08, B = 1.99 ± 0.05) and 
the scatter around the calibration curve reduces to 24% (lower 
right panel of Fig. 5). 

To obtain a goodness-of-fit estimator, the measured distri¬ 
bution is compared to the expected distribution which is com¬ 
puted from the likelihood function, i.e., from the probability 
model that describes the fluctuations. The comparison yields 
a reduced chi-square value of \ 2 /ndf = 13.8/12 for the full 
data set and x 2 /ndf = 8.43/6 for the high-quality data set. In 
particular, it shows that the estimated uncertainties of the en¬ 
ergy estimator in Sec. IV B are compatible with the observed 
scatter around the calibration curve. 


A. Uncertainties of the reconstructed cosmic-ray energy with 
the radio detector 

To determine the energy resolution of the radio detector, 
the known resolution of the surface detector needs to be sub¬ 


tracted from the combined scatter. The average (statistical) 
SD energy resolution for all air showers in our data set is 18%. 
To obtain an estimate of the radio-energy resolution we use a 
Monte Carlo study which takes into account the energy and 
zenith angle dependence of the SD energy resolution. The 
combined scatter is simulated for different radio-energy reso¬ 
lutions, according to the number of air showers and the energy 
and zenith distribution of the data set. We find that the energy 
resolution of the radio detector is 22% for the full data set and 
17% for the air showers where the core position could be de¬ 
termined in the radio LDF fit, when five or more radio stations 
have a significant signal. 

In the above calculation we assumed that the energy esti¬ 
mates from the SD and radio reconstruction are uncorrelated 
for a fixed energy. However, an anti-correlation is expected as 
radio emission originates from the electromagnetic part of the 
air shower whereas the SD signal is mostly due to muons re¬ 
sulting from the hadronic shower component [47] and which 
are anti-correlated shower parameters for a fixed cosmic-ray 
energy. In case of an anti-correlation, the estimated radio¬ 
energy resolution would be even smaller making the above 
values conservative estimates. 

Furthermore, we studied the effect of a possible bias in 
the SD reconstmcted energy for different primaries where the 
detector is not fully efficient (0.1 EeV - 0.3 EeV) and has a 
slightly different efficiency curve for the two extreme scenar¬ 
ios of proton and iron primaries [30]. We found that the effect 
is negligible for our data set. 

The uncertainty on the absolute scale of the energy estima¬ 
tor as discussed in Sec. IV C is calibrated out by correlating 
S'radio with Eqr. The method, however, inherits the uncertain¬ 
ties of the SD energy scale. This scale uncertainty is domi¬ 
nated by the FD scale uncertainty, which is used to calibrate 
the SD. It is 14% at energies > 10 18 eV [48] and increases to 
16% at 10 17 ’ 5 eV. 
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FIG. 5. (top) The radio-energy estimator SWdio as a function of the 
cosmic-ray energy Ecr. measured with the surface detector. A power 
law is fitted to the data using a likelihood approach which takes all 
uncertainties and detection efficiencies into account. Green filled cir¬ 
cles denote air showers where the core position has been determined 
in the radio LDF fit, i.e., all air showers with at least five stations with 
signal. Open circles denote events with less than five stations with 
signal and use the SD core position, (bottom) Relative energy reso¬ 
lution: The energy of the radio detector is obtained using the fit in 
the left-hand figure. The left histogram contains all air showers, and 
the right histogram contains the air showers with at least five stations 
with signal (green filled circles). The expected distribution is shown 
as a gray shaded area which is computed from the fitted probability 
model that describes the fluctuations. 


B. Precision and possible improvements of the energy 
reconstruction 

We have found that the instrumental noise and the envi¬ 
ronmental influences are not the dominant contributions to 
our energy resolution. Applying the method described to a 
CoREAS Monte Carlo data set [23, 49], including a represen¬ 
tative set of shower geometries as well as shower-to-shower 
fluctuations, but no instrumental or environmental uncertain¬ 
ties, a similar energy resolution is obtained for the same de¬ 


tector layout. 

The intrinsic limitation in the energy resolution due to 
shower-to-shower fluctuations of the electromagnetic part of 
the shower is predicted to be smaller than 10% [9, 20] and 
we expect that the current energy resolution can be further 
improved. Under the condition that the LDF samples the rele¬ 
vant part of the signal distribution on the ground correctly for 
all geometries, the energy estimator should only be affected 
by the shower-to-shower fluctuations in the electromagnetic 
part of the shower. The only additional geometric dependence 
is due to the fact that the air shower might not be fully devel¬ 
oped when reaching the ground, i.e., some part of the shower 
is clipped away. As the atmospheric depth increases with the 
secant of the zenith angle, clipping mostly affects high-energy 
vertical showers. Hence, we expect an additional dependence 
on the zenith angle. In the future, with larger statistics, this 
effect will be parametrized from data and will further improve 
the energy resolution. Also, a better understanding of the de¬ 
tector and the environmental effects, such as temperature de¬ 
pendencies, will help to improve the energy reconstruction. 

Combined measurements, such as they are possible at the 
Pierre Auger Observatory, hold great potential for future im¬ 
provements of the energy resolution due to the anti-correlation 
of the energy reconstructed with the radio and surface detec¬ 
tors. 


C. The energy content of extensive air showers in the radio 
frequency range of 30 to 80 MHz 


So far, the energy content of extensive air showers in the ra¬ 
dio frequency range of 30 to 80 MHz has only been measured 
at the Pierre Auger Observatory in Argentina. However, our 
findings can be generalized by the following consideration. 

To obtain a prediction that is independent of the location 
of the experiment, i.e., a universal formula to calculate the 
radiation energy from the cosmic-ray energy, the calibration 
function Eq. ( 6 ) can be normalized to the local magnetic field. 
We found that it is sufficient to correct only for the dominant 
geomagnetic part of the radio emission. This is because the 
increase of radiation energy due to the charge-excess emission 
is small, as constructive and destructive interference with the 
geomagnetic emission mostly cancel out in the integration of 
the energy densities over the shower plane, see Eq. (5). For the 
average relative charge-excess strength of 14% at AERA [15] 
the increase in radiation energy is only 2%. As most locations 
on Earth have a stronger magnetic field than the AERA site the 
effect of the charge-excess emission on the radiation energy 
will be even smaller. Within the statistical accuracy of the 
calibration function this effect can be neglected which leads 
to the universal prediction of the radiation energy 


E 30 SO MHz =(15.8 ± 0.7(stat) ± 6.7(sys)) MeV 


x sm a 


Ecr ^Earth \ * 

10 18 eV 0.24 G ) 


(7) 


where Eqr is the cosmic-ray energy, ar th denotes the lo¬ 
cal magnetic-field strength and 0.24 G is the magnetic-field 
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strength at the AERA site. The systematic uncertainty quoted 
here is the combined uncertainty of SVadio (28%) and the SD 
energy scale (16% at 10 17 5 eV). This formula will become in¬ 
valid for radio detectors at high altitudes because the amount 
of radiation energy decreases as - depending on the zenith an¬ 
gle - a significant part of the air shower is clipped away at the 
ground. 

Please note that in practice the 30 to 80 MHz band is used 
by most experiments. Due to coherence effects, the cosmic- 
ray-induced radio emission is strongest below 100 MHz. At¬ 
mospheric noise and short-wave band transmitters make mea¬ 
surements below 30 MHz unfeasible. From 85 to 110 MHz the 
FM band interferes with measurements. Furthermore, radio 
emission at frequencies well beyond 100 MHz can be detected 
only in very specific geometries (observers at the Cherenkov 
angle) [50], Hence, ground-based experiments exploit the 
frequency window from 30 to 80 MHz or measure in only 
slightly different frequency bands. 


VI. CONCLUSIONS 

The Auger Engineering Radio Array is the radio detector 
of the Pierre Auger Observatory. It is located within the low- 
energy extension of the Observatory where additional surface 
detector stations with a smaller spacing are present, which 
enables access to cosmic-ray energies down to 0.1 EeV. For 
the analysis presented here we only use the thoroughly cali¬ 
brated 24 LPDA radio stations of the first stage of AERA de¬ 
ployment, with data collected between April 2011 and March 
2013. 

At several observer positions, the energy deposit per area 
of the radio pulse of an extensive air shower is measured. Us¬ 
ing recent progress in understanding the lateral signal distri¬ 
bution of the radio signals, this distribution is described by 
an empirical function. The spatial integral of the lateral dis¬ 
tribution function gives the amount of energy that is trans¬ 
ferred from the primary cosmic ray into radio emission in the 
30 to 80 MHz frequency band of AERA during the air-shower 
development. We measure on average 15.8 MeV of radiation 
energy for a 1 EeV air shower arriving perpendicularly to a ge¬ 
omagnetic field of 0.24 G. The systematic uncertainty is 28% 
on the radiation energy and 16% on the cosmic-ray energy. 

This radiation energy - corrected for different emission 
strengths at different angles between shower axis and geomag¬ 
netic field - is used as the cosmic-ray energy estimator S ra dio- 
A comparison of .S' la d lo with the cosmic-ray energy of the sur¬ 
face detector reconstruction shows that it is consistent with 
quadratic scaling with the cosmic-ray energy S ra di 0 cx E B 
where B = 1.98 ± 0.04 as expected for coherent radio emis¬ 
sion. 

The calibration function is normalized to the strength of the 
local geomagnetic field. Hence, with the knowledge of the 
local geomagnetic field and a measurement of the radiation 
energy (in the AERA frequency range) the calibration func¬ 
tion can be used at any location to calculate the cosmic-ray 
energy. 

Investigating the scatter around the calibration curve and 


subtracting the resolution of the surface detector we find that 
the energy resolution of the radio detector is 22 % for the full 
data set, and 17% for the events with more than four stations 
with signal, where the core position could be determined in 
the radio LDF fit. Given the small shower-to-shower fluctua¬ 
tions of the electromagnetic component, we expect that with 
a deeper understanding of the detector and environmental ef¬ 
fects, an even improved precision in the energy measurement 
can be achieved. 


Appendix A: LDF parameters 

Table III gives the parameters used in the LDF function of 
Eq. (4). 

TABLE III. Parameters Co - Ca of Eq. (4). Co = 16.25 m and Ca = 
0.0079 m . The zenith-angle dependent values used to predict the 
emission pattern are given for zenith angle bins up to 60°. 


zenith angle 

C 0 

Ci[ 

m 

1 

C 2 [m] 

0° ■ 

- 10° 

0.41 

-8.0: 

± i 

0.3 

21.2 

± 

0.4 

10° 

- 20° 

0.41 

-10.0 

± 

0.4 

23.1 

± 

0.4 

O 

O 

CN 

-30° 

0.41 

-12.0 

± 

0.3 

25.5 

± 

0.3 

CO 

o 

o 

- 40° 

0.41 

-20.0 

± 

0.4 

32.0 

± 

0.6 

o 

o 

- 50° 

0.46 

-25.1 

± 

0.9 

34.5 

± 

0.7 

50° 

-60° 

0.71 

-27.3 

± 

1.0 

9.8: 

£ 

1.5 


Appendix B: Likelihood function 

The likelihood function (for one pair of radio signal SVadio 
and SD cosmic-ray energy estimate E’sn) has the following 
form 

UQ ^ _ 1 £SD(E’sD,0i)eRD(^SD,0 i , < i> J :) 

( radio, SD) - N 2^ £sD (^ SDji , @i ) £rd (£) SDj1 , 0,, $ 4 ) 

Xt?RD(‘S'radio|*S , (-E , SD,z), ■■•) 

X9SD-sh(-EsD|-EsD,i, ©i) . (Bl) 

The summation is performed over all events in the selected 
data set. 3 rd (SVadio IS) •••) and g SD _ sh (E SB \E,<d) are the 
conditional probability density functions, which describe the 
probability to measure a radio signal S ra dio or energy E’sn if 
the true radio signal, energy and zenith angle are S, E and 0. 
$ denotes the azimuth angle. <?RD(S ra dio) is obtained for each 
event in a Monte Carlo simulation where all reconstructed pa¬ 
rameters that influence the radio-energy estimator are varied 
within their uncertainties. £sd(E'SD, 0 ) and £rd(-Esd, 0 , 4>) 
are the efficiencies of the surface and the radio detector. The 
radio efficiency has been determined with Monte Carlo air- 
shower simulations and a full-detector simulation and depends 
on the energy, the zenith and the azimuth angle. N is the nor¬ 
malization of the function to an integral of one. 
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